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Abstract. We study the effect of non-relativistic, thermal matter in the jets of active galaxies (AGN) on the low-frequency 
non-thermal emission and the variability thereof. In matter-dominated jets, sizable quantities of gas should exist, in particular 
in the compression zones near the collision fronts that are an implicit ingredient of Fermi-type particle acceleration scenarios. 
Non-relativistic thermal gas in AGN jets noticably contributes to the optical depth at radio to infrared frequencies, and much 
less to the emission, with an efficiency that is strongly temperature-dependent. The observable flux of low-frequency emission 
is thus modulated by the temperature evolution of the thermal gas, and it can therefore display very complicated variability. For 
a particular particle energisation scenario we calculate the temperature evolution of the thermal plasma as well as the radiation 
transport of low-frequency emission, and thus derive simulated light curves at different frequencies and their typical correlation 
properties. 
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1. Introduction 

Blazars, a sub-class of active galactic nuclei (AGN), are 
among the strongest known radiation sources in the universe. 
Observations show that they exhibit strong optical polariza- 
tion, variability on all observable timescales, flat-spectrum 
radio emission from a compact core, and many of them 
emit the bulk of their luminosity in the form o f y-rays (e.g. 
iDermer & Gehrelslll995t iMukheriee etailll997h . On account 
of these measurements it is usually argued that the y-rays are 
produced in very compact regions within a relativistically mov- 
ing system, the so-called jets, thus providing a strong Doppler 
amplification of the radiation. Often, the relativistic bulk mo- 
tion can be directly observed as apparent superluminal motion 
of individual emission regions in the jets in sequences of VLBI 
observations of their radio emission. The Lorentz factors (and 
Doppler factors) thus d erived are of the order of ten for gen- 
eral samples of AGNs JVermeulen & Cohenll994\ but may be 
high er for AGNs showin g prominent gamma-ray emission (e.g. 
iHoman et alJl2002T.l2003l) . The range of Lorentz factors preva- 
lent at the time of gamma-ray emission, which presumably oc- 
curs before the emission region becomes visible at radio fre- 
quencies, is not known, but is likely higher than ten, if the bulk 
kinetic energy of the jets i s the energy reservoir for the parti - 
cle acceleration (see, e.g., iGeorganopoulos & Kazanasl2*0 03). 
These conclusions are further supported by the absence of y— y 



pair absorption and the violation of the Elliot-Shapiro relation 
jElliot & Shapir o! 1974 in several observed AGNs. 
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The unified model for active galaxies JUrrv & Padovanil 

1995) assumes that a supermassive black hole is located in 
the center of a galaxy. The black hole supposedly accelerates 
plasma to high energies, while the surrounding host galaxy pro- 
vides a steady inflow of matter, often assumed to be in the form 
of an accretion disk, to keep the system running for a long time. 
Above the accretion disk plasma is ejec ted in jets, which may 
reach lengths of s everal hundreds kpc (Begel man et alj fl984: 
ISchlickeiseJ200l. 

Though thermal radiation can be observed from many 
AGN, it is usually associated with the accretion disk and the 
medium around it. The emission from the jets is generally as- 
sumed to be entirely non-thermal on account of the spectrum 
and the variability behaviour. Nevertheless, there may be ther- 
mal plasma in the jets, in which the dominantly radiating, en- 
ergetic particles are embedded and confined. 

In this paper we will study the effect of thermal matter in 
the jets on the low-frequency non-thermal emission and the 
variability thereof. We will do so in the framework of one 
particular model of particle acceleration, in which the ener- 
getic particles are provided by the isotropization of interstel- 
lar m atter in the downstream re gion of a relativistic collision 
front JPohl & Schlickeiseii20 00, henceforth referred to as PS), 
though most of our treatment, and hence most of the results, 
are not restricted to the particulars of this scenario. 
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As we will see, the thermal gas will not necessarily mani- 
fest itself by its emission, but rather by its absorption. Both the 
emission and the absorption properties of thermal matter de- 
pend strongly on the temperature, and thus the problem at hand 
has two aspects. First, the radiation transport of low-frequency 
emission through thermal gas must be calculated, which is a 
function of the plasma temperature and density as well as of 
the spectrum of energetic particles. We then need to follow the 
temperature evolution of the thermal gas, which depends on the 
its emission and absorptions coefficients (and hence on its tem- 
perature), on the spectrum of energetic particles, and on possi- 
ble wave damping. In this study we concentrate on continuum 
emission processes. 

The spectrum of non-thermal particles obviously varies 
with time on account of the variability in the high-energy emis- 
sion. Consequently the plasma temperature will vary, and so 
will the thermal absorption coefficient. The observable flux of 
low-frequency emission is thus modulated by the temperature 
evolution of the thermal gas, and it can therefore display very 
complicated variability behaviour, the study of which is the 
subject of this paper. 

We first give a brief introduction into the PS model for ac- 
tive galactic nuclei (Sect. [2}. After this, in Sect. [5] we discuss 
the radiation transport for a general non-thermal population of 
electrons in a warm electron-proton plasma. We present numer- 
ical absorption and emission coefficients for a generic situation 
(when this is possible) as well as for the PS model (when a 
generic treatment is not possible). In Sect. |3J we discuss the 
temperature evolution in a generic warm background plasma 
for different heating processes. The results of this analysis are 
independent of the non-thermal particle spectrum, so that they 
can easily be applied to other (non-)astrophysical situations. 

In Sect. 131 and 15 .51 we present numerical results for the PS 
model and demonstrate that only a few parameters actually in- 
fluence the qualitative behaviour of a system. Section 1531 em- 
phasizes the limits of our calculations, while Sect. [6] connects 
these limits to free parameters in this situation. 

2. The basic model 

IPohl & SchlickeiseHfcOOti) have studied the kinetic relaxation 
of particles that have traversed a parallel collision front from 
the upstream to the downstream region. This pick-up process 
occurs due to scattering off low-wavenumber plasma waves 
that the picked-up particles generate through streaming insta- 
bilities. The calculation describes the first half-cycle of the 
standard shock acceleration process for relativistic, collision- 
less flows. If the outflow plasma propagating through the jet 
is dense, the downstream region provides a target for the 
ultra-relativistic protons to interact with and produce radia- 
tion through various leptonic and hadronic emission channels. 
The predicted radiation properties resemble those of gamma- 
ray blazars and thus the model offers an attractive alterna- 
tive to shock acceleration of electrons and subsequent Inverse 
Compton scattering. 

The geometry of the system is shown in Fig.[2for two sys- 
tems of reference, the jet frame (without index) and the labo- 
ratory (host galaxy) frame, for which all quantities are indexed 



lab system 




jet system 
r 




2r 



Fig. 1. The geometry visualized for two systems of reference, 
the jet frame and the laboratory (host galaxy) frame, for which 
all quantities are indexed with an asterisk. We consider a single 
cloud of thermal proton-electron plasma with density n p , that is 
relativistic ally moving through the interstellar medium of the 
host galaxy with density n* . An observer would see radiation 
from the plasma cloud that is emitted at an angle 9 with respect 
to its direction of motion. 



with an asterisk. In this report, we also use a third reference 
frame, the observer frame, which we use to calculate the ob- 
served evolution time. We will come back to this point later in 
this section. We consider the jet to consist of individual clouds 
of thermal proton-electron plasma, that are relativistically mov- 
ing through the interstellar medium of the host galaxy. Here we 
discuss the interactions of one such cloud with temperature T, 
density n p , and as a simple approximation for its spatial extent 
and form a cylindrical disk with radius r and thickness d. We 
emphasize that we treat the swept-up matter as a simple collec- 
tio n of particles, which i s in a ccordance with the earlier model 
bv IPohl & Sch lickeiser (2000), and un like the fl uid treatment 
common to shock acceleration physics JSchlickeiserl2002l) . As 
a consequence, there is no relation between the two quantities 
tip and tii, and the Rankine-Hugoniot conditions do not apply. 

The Lorentz transformation between the two frames of ref- 
erence gives the relations n\ = Tn* and n* = Trip with a relative 
Lorentz factor of T. The swept-up particles are of high energies 
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in the instantaneous downstream (jet) frame and are quickly 
isotropized JPohl & Schlickeis'eife oOO). Momentum conserva- 
tion causes the system to decelerate, while internal processes 
produce a broad distribution of primary protons, N p (y ? ), and 
secondary electrons, N e (y e ), which are more important than 
primary electrons for the parameters of interest. Unless noted 
otherwise, all relevant expressions use cgs units, all particle 
spectra N x are taken per unit volume and all equations and for- 
mulae are taken in the instantaneous jet frame. 

The isotropization of the swept-up particles effects a mo- 
mentum transfer from the interstellar medium particles to the 
jet. Therefore, the jet decelerates and the Lorentz factor of 
the jet, T, is a function of time. It follows the equations 
JPohl & Schlickeiserl2OO0l) 



71 r BjCfflp , ,„ 

r = — t~t~. — (f — i) 



M{t) 



M(t) = V 



Xoo 
dym x (n b + yN x (y,t)) 



(1) 

(2) 



x=e,p 



where M(t) is the total mass of the compact plasma, which in 
turn depends on the differential number density of the ener- 
getic particles, N(y). Such gradual deceleration has been di- 
rectly observed in th e case of the micro-quasar XTE J 1550-564 
JCorbel et al.l2002h . 

The Lorentz factors observed in astrophysical sources vary 
greatly depending on the source class and the measured veloc- 
ity indicator. While in gamma-ray burst (GRBs) they may reach 
values of several hundreds, the Lorentz factors in AGN jets are 
usually only of the order of ten, if deduced from measurements 
of superluminal motion, but may be higher for y-ray-bright 
blazars. For TeV-blazars the radiation modeling often also re- 
quires significantly higher Lorentz (and Doppler) factors (e.g. 
iKonopelko et all200 3). However, gamma-ray emission is pre- 
sumed to occur before the emission zone becomes visible at 
radio frequencies, and the corresponding Lorentz factor should 
thus be very much higher than ten, if the bulk kinetic energy 
of the jets is supposed to be the source for particle accelera- 
tion. Additionally, the physics of jet formation, and therefore 
the Lorentz factors involved with this process are not under- 
stood, but considering the supposedly extreme conditions near 
a black hole, it is highly possible that the particles which are 
ultimately injected into the jet initially possess similarly high 
velocities. For this reason, we assume an initial Lorentz factor 
To = 300. We would like to point out that our results are only 
parametrically dependent on this choice. 

In the instantaneous jet frame the differential number den- 
sity of protons and electrons follows the evolution equation 



dN ■ d , , N{y) 
ot dy T c 



(3) 



where T c is the timescale for particle losses by escape or decay, 
y represents continuous energy losses, and N is the injection or 
production rate of new particles. As initial condition at the time 
t — 0, when the jet Lorentz factor is F = Tq, we assume that 
there are no energetic particles in the jet, N(y) = 0. 

In the default configuration the background plasma sweeps 
up protons from the interstellar matter, while electrons are 



produced in inelastic collisions in the system. The dif- 
ferential proton source d ensity in the co-moving frame is 
JPohl & SchlickeiseifeOOrit 



n = — -i- Vr 2 - 1 6(y - r) 

a 



(4) 



while the (secondary) electron production proceeds through the 
reaction chain 



p + p — > + X 

TT* -»//* + V 

/i* -> e ± + 2v, 



(5) 
(6) 
(7) 



and is calculated numerically to properly account for the 
pion multiplicity spectrum in reactio n (3). Here we use the 
Monte Carlo model DTUNUC (V2.2) jMohring & RanfJl99lt 
iRanft et alJl994HFerrari et alll996ailEngel et alll997t> which 
is based on a dual parton model (Canel la et al.l Il994f> . This 
MC model for hadron-nucleus and nucleus-nucleus interac- 
tions includes various modern aspects of high-energy physics 
and has been successfully applied to the description of hadron 
production in high-en e rgy collisions (Ferrari et al.l fl9 96b: 
I Ranft & RoesleJ Il994t iMohring et ail Il993t iRoesler et alJ 

The particle losses are diffusive escape (protons and elec- 
trons), neutron escape following p+p — > n+X reactions and 
pair annihilation (electrons and positrons), 



1.67-10" 



d 2 Tn* 



sec 



= 310 15 n p 1 (exp(- 



3cctt-« p 



yl.l ■ 10 13 cm 



8(y+ 1) Vr 2 - i 



(8) 

(9) 
(10) 



(r + 4 + y- 1 )log( r + ^y 2 -\)-p{y + 3) 



sec 



(11) 



For all relevant cases the exponential in eq. [5] (where L is the 
di stance to the plasma c loud boundary) can be approximated by 
1 JPohl & Schlickeiser 2000). We approximate the continuous 
losses by 



- tei, P = 3.75 ■ 10~ 16 « p 7 



-7ei, e = 7.5 ■ 10-° ra p 



r 



sec 



-7 m-i6 (r-i) 2 -i 

-7inei, P = 7-10 n„ — sec 



-tSync = 1.1-10 



7+1 

B 2 



y 2 B 2 sec 



-TBrems = 8 • lO^ 6 « p (7 - 7~') Sec" 1 . 



(12) 
(13) 

(14) 

(15) 
(16) 
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parameter 


default value 


average magnetic field B 


1 G 


background plasma density n p 


5 • 10 8 cnr 3 


interstellar matter density n* 

j i 


0.4 cirr 3 


observer angle #° bs 


0.1° 


initial Lorentz factor T 


300 


disk height d 


3 • 10 13 cm 


disk radius r 


10 14 cm 


initial temperature T 


10 4 K 



Table 1. Standard parameters used in this report unless noted 
otherwise. For the relation between the observer angle and the 
emission angle see Eall7landll8l 



Here the subscripts el and inel stand for elastic and inelas- 
tic processes, Sync for synchrotron radiation and Brem s for 
Bremsstrahlung. The equations[21and[0]are taken from Haug 
lll988l) . eq. 18, assuming a Coulomb logarithm of 25 and ne- 
glecting the last terms. This modi fication has been m ade to be 
° ° . . I 1 i -j 

consistent with eq. |441 while in Pohl & Schlickeiser (2000) a 
value of 20 is used. It does not af fect the overall res ults in any 
visible way. Equation[2]is from Pach olczvkl(ll970l) . averaged 
over all emission directions. The inelastic energy losses, eq. 
[TH are a result of the model used to evaluate these processes 
( Cao ella et al.l 19941). Finally , the Br emsstr ahlung energy losses^ 
lHav 



are ta ken from iHavakawal dl969h and IPohl & SchlickeiseJ 
(2000). 

The timesca les for particle losses are from 

I.Tauch & Rohrlichl J 19761) (eq. US and IPohl & Sc hlickeiser 
(2000) (eq.|8]and|9}. The diffusive escape timescale used here 
is derived for the disk geometry under the condition r » d, 
which implies that effectively the charged particles will only 
escape through the top and bottom surfaces. The most impor- 
tant point here is the assumption that the energetic particles 
start diffusing outwards near the center of the jet plasma cloud. 
Depending on the Alfven speed and the initial intensi ty o f 
scattering waves this may be questionable JVainio etall2 004). 
We also want to emphasize that these equations have been 
derived under the assumption that the plasma is fully ionized 
and that the particle number is conserved. If this is not the 
case, the above equations have to be modified for ionization, 
recombination and pair production. 

The free parameters of the model are specified in table [2 
for a quick reference, where we also show the numerical val- 
ues we have adopted as standard set of parameters in this paper. 
Whenever we modify one of these values, it is explicitely men- 
tioned in the text. We have neglected all e + <?~ annihilation pro- 
cesses, since we focus on the optical-to-infrared region where 
pair annihilation does not contribute. 

All calculations in this analysis are made in the jet rest 
frame. For a comparison with observations, the photon spec- 
tra have to be transformed into the observer's frame, which, 
depending on the cosmological redshift of the AGN in ques- 
tion, can be different from the host galaxy frame. For ease of 
exposition we will assume in the following treatment that the 
observ er's frame and the host galax y frame are identical, as was 
done in Pohl & Schlickeiser (2000). 



As a reminder we here list the relations between relevant 
quantities in the observer's frame and the jet frame. While the 
observer's frame is fixed, the jet frame is only the instantaneous 
rest frame of the jet plasma at the time considered. Since we are 
interested in the long-time evolution of a compact relativistic 
object moving along a straight line under a. fixed aspect angle 
in the observer's frame, we also need to consider that the cor- 
responding aspect angle in the instantaneous jet frame is not 
constant, for the Doppler factor is changing with time. 

cos 6 + j3r 



gobs _ 



cos = 



1 + B T cos 6 
cose° bs -B r 



1 



D = 



, ,obs 



■B r cos0° bs 
1 



jobi 



df 



.obs 



F(l -B T co&e° bs ) 
Dv 

D 3 I(v" 
dt/D 



7D) 



(17) 
(18) 

(19) 

(20) 
(21) 
(22) 



3. The radiation transport 

3.1. The radiation transport equation 

To calculate the emitted photon spectra 7 Vj o for a given line- 
of-sight, it is required to know the emission coefficient j v and 
the absorption coefficient a v . The se two quantities enter in th e 
radiation transport equation (e. g . |R vbicki & Li ghtmanl 1 9 8H) . 

^=S V -I vfi (23) 
dr 

via the source function S v — jil 2/ a i an d tne optical depth 

a dx . (24) 



The well-known general solution to this equation are for a path 
length of s 



/v,o(r) 



= Jv 



[1 
■ s 



■ exp(-r)] 



for r <K 1 
for t » 1 . 



(25) 



In our calculations we have assumed that the radiation coeffi- 
cients j Y and a v are independent of location, i.e. the supposedly 
cylindrical cloud of plasma in the jet is homogeneous, so that 
the optical depth reduces to 



t s = a ■ s 



(26) 



These expressions are the solution of the radiation trans- 
port equation for a single line-of-sight through the jet plasma. 
We still have to integrate over the entire emitting surface, as 
well as respecting the exact geometry of the system, since the 
photon path length is not constant over the apparent surface of 
the emission region. The total spectral power is then 



dA(Q) 7 V)0 (fl) 



(27) 



with the apparent surface element dA and the emission direc- 
tion O. 
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Since a real cloud of plasma will have a complicated (and 
generally unk nown) geometry rather tha n the simple shape 
considered bv lPohl & Schlickeiserl l Eo OO) and in this report, it 
may be sufficient to arrive at an approximate solution to Eq|27l 

Let us start with the most simple approximation: we assume 
that the path length, s, is independent of the emission point both 
for the front surface, A t , and the side surface, A s . Then for the 
front surface 



Sf - d/cosO, 
and 

P v ,f(£l) A f cos 9 I v fi(d6) 



nr cosfl5\ 



1 - exp 



ad 
cos 9 



(28) 

(29) 
(30) 



On account of the assumed thin-disk geometry, i.e. d <K r, the 
contribution from the side surface is significant only for emis- 
sion angles 9 w tt/2. The choice of a thin disk for our system 
introduces an additional problem, for the photon path length 
between two opposite points on the side surface is not constant; 
we use s =2 r and obtain 



1e-05 



1e-10 - 



5 1e-15 



1e-20 



1e-25 



synchrotron 
free-free 
synchrotron (no Razin) 



-5.5 



-4.5 -4 -3.5 
log 10 E v [eV] 



-2.5 



Fig. 2. The absorption coefficients a for the free-free pro- 
cess and the synchrotron process for the jet model of 
IPohl & SchlickeiseH fcOOOh with standard parameters (see ta- 
ble[Q. We display the situation after one hour in the observer's 
frame, at which time the plasma temperature is T« 3 • 10 4 K. 
Apparently the system is dominated by free-free absorption. 



s s ~ rj sin#, 
and 



P v ,s(&) = A s sm9I Yfi (d&) 



2d r sinflSv 



(31) 

(32) 
(33) 



which turns out to be a fair approximation to the exact solutions 
of the problem, which are derived in appendix lAl 

A better approximation is afforded by a box geometry in- 
stead of the original cylinder. The side lengths would be a — 
\fiir to conserve the volume. If one assumes that the line-of- 
sight to the observer is in the plane of two of the side sur- 
faces, EqEJcan be solved exactly. The resulting expressions 
for the observable photon spectra are presented in appendixlBl 
The box geometry turns out an excellent approximation to the 
exact solution for the initially assumed disk geometry, with a 
maximum error of 3% at 9 — 60°, while the simple constant- s 
approach introduces an error of 20-30%. 

3.2. The emission coefficients 

In this report we consider two different contributions to 
the radiation processes, synchrotron emission and thermal 
bremsstrahlung (free-free emission). The bremsstrahlung co- 
effici ents for a qu asineutral thermal hydrogen plasma are 
JRvbicki & Lightmanl 19851) 



yln 1 acrj he mc 2 (2kT 
jbr = — \l -r=- In 



T 3/2 



kT 



hv 



hv 

exp|-- 



(34) 



a hr = 



T 2 3 

2n a o~t c 
(27r) 5 / 2 v 3 



mc 



In 



2kT 
hv 



1 - exp - 



hv 
kT 



(35) 





1e-18 r 




1e-19 r 




1e-20 r 


tfl 




N 




X 

CO 


1e-21 r 


E 









■S3 


1e-22 r 












1e-23 : 




1e-24 r 




1e-25 




synchrotron 
free-free 
synchrotron (no Razin) 

-4.5 -4 -3.5 
log 10 E v [eV] 

Fig. 3. The emission coefficients j for the free-free process and 
the synchrotron process after one hour observed time for stan- 
dard parameters, as in Fig|2] The free-free emission is only vis- 
ible because of the Razin effect. For « p = 5 • 10 s cm~ 3 the 
plasma frequency corresponds to E v » 8 • 10~ 7 eV, so we cover 
the entire valid energy spectrum. 



The synchrotron radiation coefficients for an ensemble of 
electrons with differential number density N{y) can be calcu- 
lated as JRvbicki & Lightmanl 19851) 



1 r 

hy = J dyN(y)p sy (v) 



(36) 



a - = -^r dy ^- i)p - iv) U0i (3?) 



with the spectral power per electron 
V3e 3 fi 



Psy(v) 



mc- 



F\ V - 



(38) 
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Fig. 4. Emitted photon spectra after one hour of observed time 
in the jet system for the standard parameters (6° bs * 0.1°, see 
table Q and T ~ 3 ■ 10 4 K. The thermal Bremsstrahlung dom- 
inates at low energies only because of the Razin-effect, which 
suppresses the synchrotron spectrum here. 



where 



3eB 
Anmc 



y 



F(x) * 1.8 x 3 exp(-x) 



(39) 



In a dense plasma with plasma frequency v p the spectral 
power is modified at low frequencies by the Razin e ffect 
ilRvbicki & Lightmanlfl98l ICnisius & SchlickeiseJfl988i) on 
account of the modified dispersion relation for electromagnetic 
waves. Then 



Psy,Razin(v) 



y (vR i/2 ) 



where 



R = 1 



+ m 2 



9- 10 J 



n v Hz. 



(40) 



(41) 



The Razin effect creates a cutoff in the photon spectrum at low 
energies, whereas at high energies the emission spectrum is un- 
changed. 

The radiation coefficients for the standard para meters (see 
table □ in the jet model of IPohl & Schlickeiserl l l2000l) are 
displayed in Fig[2]and Fig[3] The synchrotron coefficients are 
based on the numerically derived electron spectrum in the sys- 
tem after one hour of observed time, while the free-free radia- 
tion coefficients depend only on the temperature, which is ap- 
proximately 3 • 10 4 K at this time. The kinematic of the pion and 
muon decays is unfavorable for the production of secondary 
electrons with an energy less than approximately 100 MeV. 
Thus, essentially only highly relativistic electron are generated. 
The energy loss rate below 300 MeV is dominated by Coulomb 
and ionization interactions, i.e. is independent of energy, so that 
in a steady-state spectrum the cooling tail towards lower ener- 
gies would be flat. The synchrotron spectrum in Fig.[3]is corre- 
spondingly inverted. 

The non-thermal emission spectra do not significantly 
change when we modify the parameters, except for an approx- 
imately linear dependence on the density of the interstellar 



medium, so the results displayed here are a good estimate in 
most situations. The free-free absorption, and to a smaller ex- 
tent the emission, is strongly reduced in efficiency for higher 
temperatures. As we shall see later, the plasma temperature 
can rise from T — 3 • 10 4 K, which nearly is a lower limit, 
to 10 9 K or more. Modifying the temperature from 3 ■ 10 4 K 
to this value reduces the absorption coefficient by 7 orders of 
magnitude, which makes the free-free absorption less efficient 
than the synchrotron absorption. The temperature, for which 
both are approximately equal, is around 10° K. 

It thus appears that the optical thickness is mostly domi- 
nated by free-free absorption, while the emission is dominated 
by either synchrotron radiation or bremsstrahlung, depending 
on the frequency. It is important to include the Razin effect in 
our calculations, for, without the Razin effect, the synchrotron 
emission coefficient would not have a cut-off at low frequen- 
cies, and free-free emission would be completely negligible. 

An interesting possibility is that a synchrotron maser may 
operate in the jet plasma. Maser activity requires a negative 
absorption coefficient (a < 0), indicating that stimulated emis- 
sion is more important than the spontaneous absorption. If the 
optical depth r <SC — 1, the emitted photon intensity can be very 

high; 

I 1 j -i 

Crusius & Schlickeiser ( 1988) have shown that for the syn- 
chrotron process this is only possible when the electron spec- 
trum has a slope of at least 2, i.e. the distribution function is 
inverted. The numerical calculations of the inelastic proton de- 
cay chain seem to indicate that this is not the case, at least 
in the energy regions covered in this report. The condition for 
synchrotron masering is a < 0, and since the evaluation of a 
requires an integration over the entire radiating particle distri- 
bution (ea. l37t . this condition must be valid for a large energy 
interval, and a locally hard spectrum will not be sufficient for a 
maser to operate. 

4. Simple photon spectra 

In this section we discuss radiation spectra for a constant tem- 
perature with a view to understand the characteristic properties 
of the escaping low-frequency emission. After one hour of ob- 
served time, the temperature is T ~ 3 ■ 10 4 K. 

Photon spectra calculated on the basis of the exact disk ge- 
ometry are presented in Fig.@] The spectrum can be separated 
in an optically thick region and an optically thin region. The op- 
tical depth itself is dominated by the free-free absorption pro- 
cess. The emission process is dominated by the synchrotron 
process, for which the Razin effect causes a drop off below a 
certain photon energy, so that at very low frequencies the free- 
free emission becomes visible. 

The spectra shown here are calculated in the jet frame, and 
thus they need to be transformed in the observer system for a 
comparison with data. For blazars, this would result in a fre- 
quency shift by about two orders of magnitude, so the turn- 
over frequency, at which the emission transitions from opti- 
cally thick to optically thin conditions, would be located in the 
near-infrared. 

Next, we calculate the emitted photon spectra for different 
observation angles f5 obs and show typical results in Fig|5] For 



M. Siewert et al.: The impact of thermal gas in AGN jets on the low-frequency emission 




1e+14- 



1e+12- 



6=5° 

6 o =0.008° 
r=264 




~\ 1 1 1 1 r 

1 h obs. time 

10 h obs. time 
100 h obs. time 



0=177° 
e°=7.28° 
r=299 (297 after 100 h) 



-6 -5.5 -5 -4.5 -4 -3.5 -3 -2.5 -6 -5.5 -5 -4.5 -4 -3.5 -3 -2.5 -6 -5.5 -5 -4.5 -4 -3.5 -3 -2.5 -2 

E[eV] 

Fig. 5. Emitted photon spectra in the jet frame after one hour of observed time for different emission angles, based on an exact 
treatment of the geometry. For high fl° bs , the slow evolution of the synchrotron part of the spectrum has been reproduced. 



very small aspect angles 6 10 < 1°, the jet system appears to 
rapidly evolve; in fact F « 270 after one hour of observed time. 
This value is still very high compared to typical Lorentz factors 
of 10, but one should note that it takes months for an outburst to 
become visible at radio frequencies, where most measurement 
of jet velocity are made. 

The characteristic aspect of these spectra is the presence of 
both the thermal and the synchrotron components of the spec- 
trum. Variations on the high-frequency part of the spectrum, 
most notably a shift of the turn-over frequency, are visible, 
which result from the different aspect in the jet frame, i.e. the 
different optical depth for a frontal and side view of the emis- 
sion region in the jet frame. 

For larger observation angles (6" bs > 1°), the system 
evolves much slower, so that T « Fq after a single hour. 



5. The evolution of the temperature 

Throughout this paper we have assumed that LTE is valid, that 
the jet plasma is always non-relativistic, and that the electron 
and proton components have the same temperature, i.e. the 
internal equilibration processes are faster than the energy ex- 
change with the environment. 



5.1. The cooling process 

As we have seen, the temperature of the jet plasma is a de- 
cisive parameter for the radiation yield at low frequencies. 
The plasma is subjected to a number of heating and cooling 
processes, and therefore the temperature will not remain con- 
stant. We use a balance equation to follow the variations of the 
plasma temperature, 



t = -V(y~l)(e~A), 



(42) 



with the adiabatic index y (=5/3 for a non-relativistic hydro- 
gen gas), the heating rate, e, and the cooling rate, A. For the 
cooling rate we use an analytical approximation (see fig. |6} of 
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Fig. 6. The approximated cooling function according to Eal43l 



the standard interstellar coo ling function (Dalg arno & McCravl 
ll972URavmond et alll976l) . 
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for 10 5 < T < 4 • 10 7 

for T > 4 ■ 10 7 



(43) 



where the temperature, T, is in units of Kelvin. The error in- 
troduced by this approximation is not significant, for other un- 
certainities such as the real geometry have a much greater im- 
pact on the final results. A possible exception is the jump at 
T = 10 4 K that results from electron-impact excitation and 
ionization of H i, but may be modified here on account of the 
collisional excitation and ionization by the energetic particles. 

5.2. The Coulomb heating process 

The heating rate is more difficult to calculate, since there 
are several concurring contributions. For a completely ion- 
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ized plasma, the heating rate from elastic Coulomb scatter- 
ing of a distribution of suprathermal electrons or protons, 
N(y), in a thermal plasma of tempera ture T is given by 
(Butl er & Buckingham! 1961 lHau J 19881) 



r°° 3 

dy N(y) — cr T m e c 3 n p In A 
Ji l P 



(44) 



where In A « 25 is the Coulomb logarithm. This expression 
is valid for protons and electrons (this follows from Eg 1 131 and 
I12L so that usually the non-thermal particles with the greatest 
number density dominate this heating process. In most practi- 
cal situations in our model, this is the secondary-electron pop- 
ulation. 

If only relativistic particles are considered, y » 1 and /3 - 
1; then the heating rate scales linearly with the total number 
density of high energy particles, N tot . 



e =a - cry- m e c n p Mot In A 



5.3. The synchrotron heating process 



(45) 



Heating is also effected by absorption of non-thermal emission. 
In our case, the dominant process is free-free absorption of 
synchrotron photons. The synchrotron absorption coefficient is 
usually small compared to the free-free absorption coefficient, 
and hence the emitted radiation is almost exclusively absorbed 
by the free-free process. The heating rate can be calculated as 
an integral over the deficit in total luminosity caused by the 
absorption. 



(46) 



The spectral power, P V (Q), is given by Eq|^]and strongly de- 
pendent on the emission angle, thus necessitating the integral 
over 0. In appendix|A]we will discuss ways to arrive at approx- 
imate solutions to P v (£2) with various degrees of accuracy. 

5.4. Comparing both heating processes 

It may be useful to identify parameters, for which Coulomb 
heating is always more efficient than heating by absorption, in 
which case solving Eq|46|would not be necessary. For that pur- 
pose, we neglect the Razin effect, for the very low-frequency 
photons do not significantly contribute to the total heating rate. 
We also note that the optical depth is strongly dependent on 
frequency, thus effectively dividing the spectrum at the turn- 
over frequency v+, for which r = 1, in a low-frequency part 
at v < v+, in which essentially all photons are absorbed, and a 
high-frequency part at v > v+, in which essentially all photons 
escape. 

We have noted before that free-free absorption is generally 
more efficient than synchrotron absorption. The turn-over fre- 
quency v+ is determined by the absorption coefficient (Egl35t 
and the path length s. Setting the logarithmic term to the con- 
stant value L we find in the limit hv <k kT 



10 14 cm 



10\10 6 K/ 



(47) 



With these approximations eq. |46| reduces to a much simpler 
form. Using eq. [36] and changing the order of integration we 
obtain 



e < I dv A-n j y 
Jo 



(48) 



dyN(y) dv Psy (v)= dy N(y) f(y) SY . (49) 
Ji Jo Ji 

This integral has the same mathematical form as that in Eg 1441 
and Eq|^] implying that is may be sufficient to compare the 
weight functions for the heating processes, /. 

For Coulomb heating by relativistic particles, the weight 
function, / e i. e , follows directly from eg. 1451 



/(y)el,e -6-10 



-19 . 



i p erg cm 3 s 1 



(50) 



Since electrons and protons have the same weight functions, 
/(y)ei,e = f(y)et,p> this applies to protons as well, which are 
usually less numerous that electrons and therefore won't dom- 
inate the overall heating process. 

The synchrotron weight function is more difficult, for the 
integral over v in Eg|49l 



/(r) 



- r 

sy Jo 



dv p sy (v) 



= 1. 



V3V 



'I 



dv 



U/3 

— exp( ) 

V c j V c 



(51) 



requires a numerical treatment for general v + . However, an ap- 
proximate solution can be derived. For that purpose we com- 
pare the electron energy loss rates for synchrotron emission 
(Egll5t and Coulomb scattering (Egll3>. For electron Lorentz 
factors, y, smaller than the critical value 

y_ = maxCl.^.B)) (52) 

g (n p , B) = ^682-Jm^ * 0.024 (|) ' , (53) 

the electrons lose more energy by Coulomb scattering than by 
synchrotron radiation. This implies that Coulomb interactions 
are the dominant heating process, whatever the optical depth of 
the jet plasma. Only at Lorentz factors y > y_ will the elec- 
trons dump a significant fraction of their energy into the radi- 
ation field. At these energies, the electrons will mainly radiate 
at frequencies v > v_, where 



(54) 



i.e. under optically thin conditions for the standard parameters 
(see table[T}. So 



» v+ 
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which means that only a small fraction of the synchrotron emis- 
sion will be absorbed by the plasma, and thus Coulomb heat- 
ing should always be more efficient than heating by absorp- 
tion. A careful evaluation of Eq|^ using a Taylor expansion 
in v + «: v c indicates that the numerical factor on the RHS of 
EqI55I should be closer to unity, but, on the other hand, that 
equation does not consider the heating by energetic protons. 
We may therefore conclude that synchrotron heating can be ne- 
glected, if the condition (1551 is met. 

5.5. The balance of heating and cooling 

The evolution of the temperature is closely tied to the evolu- 
tion of the heating and cooling rates. For a homogeneous in- 
terstellar medium the system will tend to attain a quasi-steady 
state. Consequently, the heating rates rise with time until they 
reach an asymptotic value. The temperature then evolves until 
the cooling balances the heating. 

The cooling rate is a strong function of the temperature, 
with a huge jump at 10 4 K as the dominating feature at low tem- 
peratures (see fig|6j. To reach higher temperatures, the heating 
rate must be large enough to compensate the hydrogen cool- 
ing, or the system will stay at 10 4 K, where the ionization frac- 
tion of hydrogen may be low, which results in an even lower 
turn-over frequency, so that the entire emission may be opti- 
cal thin. Additionally, it is likely that under the conditions pre- 
vailent in the jet plasma, the temperature dependence of the 
ionization fraction is weake r tha n in the LTE ca l culatio ns of 
Dal garno & McCravl i 19721) and IRavmond etaD dl976l) . thus 
causing a slower rise of the cooling curve. Nevertheless, the 
system will have a stable region around 10 4 K for a wide range 
of heating rates. 

When the heating rate becomes larger than the maximum 
of the cooling function at 10 5 K, the temperature will increase 
to more than 10 8 K, beyond which bremsstrahlung cooling may 
again balance the heating. Extrapolating the cooling function to 
temperature T > 10 8 K is problematic, though, for one quickly 
arrives at relativistic temperatures, kT > m & c 2 , at which most 
of our assumptions and formulae are no longer valid. At these 
temperatures, it is especially required to include pair produc- 
tion into our equations, which increases the particle number 
in the background plasma while reducing the amount of energy 
available for the individual particle, thus efficiently limiting the 
temperature of the system. In addition, the concept of a stable 
isotropic Maxwell-Boltzmann distribution is not well justified 
in these regions (Stepne \ll982l) . 

For the standard parameters (see table this overheat- 
ing does not happen (see fig. [7}, while for n* = 0.8 crrT 3 the 
heating processes are too strong to keep the temperature non- 
relativistic. The Fig. also demonstrates that the synchrotron 
contribution to the heating rate is quite small; numerical eval- 
uations show that this process contributes about 10 % of the 
Coulomb heating rate, which is in good agreement with our 
estimate in the preceding section. Finally, the asymptotic cool- 
ing is a direct result of the deceleration of the system, which 
decreases the number of swept up particles (eq. 0J. At t — 0, 
when there are no non-thermal particles present, the heating 
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Fig. 7. The temperature as a function of time for the standard 
parameters as well as n* = 0.8 cm 3 , for which the temperature 
becomes larger than 10 9 K, where our model is invalid. 



rates equal zero, and the system usually cools down to the sta- 
ble point at 10 4 K. It is possible to prevent this with a careful 
choice of initial parameters, which results in stable tempera- 
tures of T = 10 8 K. 



6. Results 

6.1. Modifying the initial parameters 

In the last sections we have seen that the temperature of the 
background plasma can be used to estimate the stability of our 
model. For this reason, we have made a parameter study to 
approximate the regions in which we get reliable results. We 
have discovered four different regions: cold (T < 10 4 K), warm 
(10 4 K< T < 10 5 K), hot (10 8 K< T < 10 9 K) and forbidden 
(T > 10 9 K). The region between 10 5 K and 10 s K is unstable, 
which is a result of the specific form of the cooling function. 
The region between 10 8 and 10 9 K is stable only for very spe- 
cific initial parameters, as outlined in the last section. 

The particle number densities « p and n* strongly dominate 
the system; the different regions in the n p -n*-plane for the stan- 
dard parameters (see tabled are displayed in Fig[8] Another 
parameter which determines the temperature evolution is the 
size of the system. It turns out that on account of the parti- 
cle number sweep-up rate (eq|4j, which is proportional to the 
sweep-up surface divided by the total volume of the system, 
only one of the size parameters really contributes to this point, 
which can be identified as the average optical path for a generic 
system. For the specific disk geometry, this parameter is iden- 
tical to the thickness d; then for the standard parameters the 
boundary points are 



d co , 



Id-warm 



n—hot 



2.1 ■ 10 13 cm 
4.1 ■ 10 13 cm 



(56) 
(57) 



Since the secondary electrons dominate the heating process 
(see Sect. 15. 21 . and their production rate can not be treated an- 
alytically, these results can not be derived from Fig.|8] 
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Fig. 8. The three qualitatively different temperature regions as 
a function of the particle densities and the standard parameters 
(To = 10 4 K). The dotted lines are for T = 600. 

The initial Lorentz factor T(f = 0) modifies the interstellar 
matter density in the jet frame, n\ - F(t)n*. This in turn mod- 
ifies the number of particles in the system and therefore the 
heating rate by elastic scattering. However, it does not modify 
the form of the boundaries of temperature regions in the n p -n*- 
plane, but only their location. To demonstrate this, in fig.[8]we 
also display the temperature regions for To = 600. 

The observer angle #° bs only modifies the observed time re- 
quired for the system to evolve, so the asymptotic behaviour 
remains unchanged. The magnetic field strength B has prac- 
tically no influence on the temperature, since the synchrotron 
heating process is negligible. 

As mentioned earlier, by modifying the initial temperature 
To together with the interstellar matter density n*, we are able 
to establish stable 'hot' temperatures. In most cases, however, 
the initial temperature will not modify the system. 

From this we conclude that our model is stable only on a 
rather small region in the initial parameter space. In this case, 
there are three asymptotically stable temperature regions, and 
it may be possible to fix some of these parameters by obser- 
vations. In all other cases, the initial assumption of a non- 
relativistic thermal component in the jets breaks down, al- 
though the sweep-up mechanism itself remains valid. 

6.2. Time-dependent variations of free parameters 

In a real situation the extragalactical matter density n* is not 
constant over long timescales and distances. The plasma in the 
jet may cross dense clouds of hydrogen or enter a local fluctua- 
tion. For this reason we have investigated the temperature evo- 
lution under some simple modifications of our sweep-up rate. 

First, we have considered a simple step in the ambient mat- 
ter density, 



N dd = #(l+fl0(f!-f)). 



(58) 



10 K, although this situation will still be quite rare. Generally, 
this modification simply produces a break in the temperature 
evolution. 

Secondly, we have investigated the effect of local density 
fluctuations. We have studied a situation in which the system 
periodically passes through a dense cloud of material, that is 
embedded in low-density gas, where for a single period t the 
sweep-up rate is 



Npeokif) = 



a for t < c 
b for c < t <t 



(59) 



Such a situation may possibly be found when the jet leaves the 
immediate host galaxy and enter the intergalactic medium. In 
this case it is slightly easier to produce stable temperatures at 



with a » b and c the fraction of the period where the heat- 
ing is strong. We have normalized these parameters in a way 
that the interstellar matter density averaged over a period is 
n*. Under this periodically peaked scheme, it is again easier 
to keep the system stable at 10 8 K without causing overheating 
as described in chapter l5~5l 

6.3. Light curves and correlations 

We have then calculated light curves for several different en- 
ergies, which we present in Fig.|5]for the periodically peaked 
injection and the standard parameters. Unlike the other photon 
spectra presented in this report, these curves are taken in the 
observer frame, in which, on account of the deceleration of the 
jet, the frequency corresponding to a fixed jet-frame frequency 
varies with time, although the modification usually is small. 

We see that the variability is visible, but that the average 
behaviour of the system remains unchanged. Additionally, the 
variability only affects the optical thick region, which is a result 
of temperature variations. In the optical thick region, we have 
/ oc j/a (ea.l25> and a « athermai, which in turn causes F v to de- 
pend on the temperature. Since the heating rate is proportional 
to the total number of energetic particles, a sudden increase in 
this number causes a similar modification of the temperature of 
the jet plasma, which again modifies the optical thick emission. 
Because of j w j sy , the optical thin emission (see ea. !25t does 
not demonstrate a similar behaviour. For this reason our model 
is able to reproduce different variability in different freqnency 
bands. 

Next, we have calculated correlation coefficients for this in- 
jection profile, using light curves in the observer frame between 
the infrared, optical and ultraviolet, from the plasma frequency 
up to a few 10 eV, with small and large distances between the 
considered frequencies. 

It turns out that the lightcurves at different frequencies are 
almost perfectly correlated (> 0.9), even if we compare free- 
free-dominated with synchrotron-dominated parts of the spec- 
trum. To understand this behaviour, we first note that the elec- 
tron distribution is dominated by the diffusive escape process 
(eq. |8j, while the continuous losses are too inefficient. For 
high-energy electrons the particle spectrum is exclusively de- 
termined by the Lorentz factor of the jet, F. For this reason 
all light curves produced by the electrons are highly correlated 
with each other, even between the free-free and synchrotron 
processes, since the free-free emission depends on the temper- 
ature, which again is determined by the high-energy electrons, 
similar to the synchrotron emission. Given that the optically 
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Fig. 9. Light curves for three different energies in the observer 
frame. Here we have adopted the periodic peak injection in- 
troduced in eq. [59] with the parameters a/b = 50, a period of 
r = 1000 s observed time and c = 100 s, i. e. 10% of the pe- 
riod time. Significant variations occur exclusively in the optical 
thick regime; the redshift of the source has been set to z — 0.5. 



thin synchrotron emission shows only a very weak response to 
the variations in the injection rate, most of the correlation sig- 
nal is presumably caused by the common secular trend in all 
light curves that is caused by the deceleration of the jet. 



7. Possible corrections 

In this section we briefly discuss possible further aspects of the 
problems investigated in this paper. 

7.1. The electrostatic instability 

Recently, IPohl et al. (2002) have investigated the electrostatic 
instability in a relativistic beam of electrons and protons and 
its impact on our model. Their conclusion was that the sweep- 
up spectrum effectively gets modified from a peak structure at 
y — T to a plateau distribution up to y — F, while keeping the 
number of collected particles and the low-energy photon spec- 
trum unchanged. The energy lost in this process may contribute 
to the heating of the background plasma, but there is no reliable 
estimate available of the exact amount of energy transferred. 
Observations hint that only a tiny fraction of the total energy 
lost by this mechanism can ever contribute to the heating. 

Since this process in principle involves a huge amount of 
energy, it might dramatically change the evolution of the tem- 
perature. However. IPohl et all J2002I) have only considered the 
initial situation of a cold electron-proton beam penetrating cold 
thermal plasma. The presence of isotropic energetic particles 
should strongly reduce the growth rate of electrostatic turbu- 
lence, so that at least in a quasi-steady state heating by damp- 
ing of electrostatic waves should be much less severe than sug- 
gested by the asymptotic energy loss of the incoming particle 
beam. 



7.2. The initial conditions 

Currently, our initial conditions assume that there are no non- 
thermal particles present. Although the mechanism responsi- 
ble for cr eation and stability of the jets is not complete ly un- 
derstood JUrrv & Padovanill 19951 iBegelman et alj fl984). it is 
likely that there is no time zero, at which a plasma cloud is 
relativistically expelled without containing energetic particles. 
Because of this, the start-up phase visible in our results would 
be a consequence of our particular treatment of the problem, 
and not a physical phenomenon. 

In section |6] we have found that only a few of our parame- 
ters significantly modify the asymptotic behaviour. The average 
system size only affects the observed photon spectra around the 
turn-over frequency, where r = 1, and the particulars of the ge- 
ometry are unknown anyway. The initial particle distributions 
won't modify the asymptotic temperature evolution either, be- 
cause the presence of a moderate amount of non-thermal parti- 
cles only modifies the time required to reach a balance of gains 
and losses. Variations of the matter densities have already been 
investigated in section [6] All other initial parameters do not 
modify the general behaviour of the system. 



8. Summary and conclusions 

In this paper we have investigated the evolution of thermal 
plasma in AGN jets and its impact on the optical-to-infrared 
photon emission. In matter-dominated jets sizable quantities 
of gas should exist, in particular in the compression zones 
near the collision fronts that are an implicite ingredient of 
Fermi-type particle acceleration scenarios. We conduct our 
stud y in the fr amework of the channeled outflow model of 
( Pohl & Schlickeisei 2000), who have studied the kinetic relax- 
ation of particles that have traversed a parallel collision front 
from the upstream to the downstream region. This pick-up 
process occurs due to scattering off low-wavenumber plasma 
waves that the picked-up particles generate themselves through 
streaming instabilities. The calculations, thus, describe the first 
half-cycle of the standard shock acceleration process for rela- 
tivistic, collisionless flows. If the outflow plasma propagating 
through the jet is dense, the downstream region provides a tar- 
get for the ultra-relativistic protons to interact with and produce 
radiation through various leptonic and hadronic emission chan- 
nels. 

Non-relativistic thermal gas in AGN jets no tic ably con- 
tributes to the optical depth at radio to infrared frequencies, and 
much less to the emission, with an efficiency that is strongly 
temperature-dependent. Assuming that this plasma is in a ther- 
mal equilibrum, we have calculated the temperature evolution 
resulting from the competition of radiative cooling and heat- 
ing by Coulomb processes and absorption of non-thermal emis- 
sion. 

Similar to the well-known results for the structure of the 
interstellar medium in Galaxies, we find that the stable regimes 
exist for temperatures between 10 4 K and 10 5 K, and around 
10 8 K. Below 10 4 K the ionization fraction will be small, and 
the optical depth is modified. Above 10 8 K, the thermal parti- 
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cles reach relativistic velocities, for which our model begins to 

breakdown. 

In the model of IPohl & Schlickeiserl l|2000), short-time 
variability at low energies arises on account of density fluctu- 
ations in the upstream medium. Consequently the plasma tem- 
perature will vary, and so will the thermal absorption coeffi- 
cient. The observable flux of low-frequency emission is thus 
modulated by the temperature evolution of the thermal gas, and 
it can therefore display very complicated variability behaviour. 
For simple density profiles of the interstellar gas in AGN host 
galaxies, we have calculated the temperature response of the 
thermal gas in the jet, and have then derived light curves at dif- 
ferent frequencies. For sufficiently long observing times, all of 
these light curves turn out to be strongly correlated with each 
other, independent of the proton injection scheme, only on ac- 
count of the deceleration of the jet. 
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Appendix A: Solving the radiation transport 
equation with the exact cylinder geometry 

The most general solution to the radiation transport equation 
JZDis 



dA(l - exp(-r(*, Q)) 



(A.l) 



where the optical thickness depends on the emission point x 
and the emission direction SI. The optical thickness needs to be 
calculated with respect to the geometry of the system. In the 
case of the cylinder geometry this expression can be evaluated, 
with several one-dimensional integrals remaining which must 
be treated numerically. 

To solve equation IA.1I we treat the emitting surface in 
Cartesian coordinates. (If we used polar coordinates, we would 
end up with two-dimensional numerical integrals.) We define 
our coordinate system as presented in fig lA.ll Then the prob- 
lem reduces to an integral over all impact parameters a and an 
integral over all 'slices' x = -w(a)..w(a). We will evaluate the 
integrals over the 'slices' first. 

It is useful to define the following quantities, related to the 
selected coordinate system (see fig. IA. 1 1 , which follow from 
elemental trigonometry. 



2r 

- arctan — 
d 



a c (0) = yjr 2 -d 2 1 At&n 2 
w(a) = 2 Vr 2 — a 2 
s(0) = dtan6 



(A.2) 

(A.3) 
(A.4) 
(A.5) 



The angle C is the critical emission angle, above which some 
contributions to the integral vanish. a c (ff) is the critical impact 
parameter where the mathematical form of the slices changes 
(see fig. lA.2l for a sketch of this). Above = c c this expression 
no longer makes any sense. w(a) is the 'width' of the system for 
an impact parameter a and determines the integration limits for 




a c(9) 



emission direction 



Fig. A.l. Definitons used in the text to calculate the integral 
over the emitting surface. Top view. 

the variable x. Finally s(ff) is used to divide the slice integral in 
the two different regions called / and // in fig. IA.2l With these 
conventions the integral splits in 



lv - h+n + hu+iv 

= 2( f 1 Si(6)da+ f 

J() Ja c 



S 2 (0)da) 



(A.6) 
(A.7) 



where S; are the two types of slices that are possible in the 
system (see fig. lA.2> . 

The first type of slice is the inner part of the disk (for small 
a), where some of the emitted photons (I) see a constant path 
length, where only the border region (II) is modified. Then the 
integral over S \ can be solved exactly, and we get 



S i = (w(a) - s(6)) cos 0(1 - exp(-rj/ cos 9)) 

^j(ff) cos 9 j 



+2 



f 

Jo 



( 1 - exp(a 



cos 9 s(0) cos 



))dx 



— (1 - exp(-Trf/ cos 0)) ■ 



sin cos 

((w(a) - s(0)) cos - 2 ) 

a 

+2s(0) cos 0. 



(A.8) 



(A.9) 



To solve the integral over x we have used the parameterization 
t(x) = T max • x/x max , w ith < x < x max (see fig lA.3l >. Inserting 
this expression in lA.7l and integrating over a gives us 

I] = 2(1 - exp(-T^/ cos 0)) ■ 

(cos 0(a c w(a c )/2 + r 2 arcsin a c /r) - a c d sin 

a c sin cos 
-2- ) 



+4a c d sin = L 



(A. 10) 



which is the exact solution of the radiation transport equation 
for the region defined by a < a c . 

The second type of slice, S2 contains expressions that can 
not be solved analytically. The general structure of the integrals 
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a < a. 



w(a) 




a > a r 




a = a. 




d 



w(a) 



w(a) 



Fig.A.2. Possible configurations of the integral over x as a 
function of the other parameters. 




photon emission 



apparent emitting 
surface (symmetry) 



s(a) cos(t) 




d/cos(theta) 



x/(s(a) cos(t)), x=0..s(a) cos(t) 



Fig. A.3. Detailed view of the geometry in the x-direction for 
small 'impact parameters' a. 




apparent emitting 
surface (symmetry) 



cos(t) , 




w(a)/i 



w(a)sin(t)), x=0..w(a)sin(t) 
photon emission 



Fig. A.4. Detailed view of the geometry in the x-direction for 
large 'impact parameters' a. 



sin 9 cos 9 
((d - h{a)) sin6» - 2 ). 



(A. 12) 



Inserting this in ea. lA.7l results in non-analytical integrals, 



h = cos 9(2 arccos a c w(a c )) 

r 



2 cos 9 

+2 sin 9(r - a c )(d ) 



+2 cos 6* I w(a)exp(- 



sin(0) 



)da) 



(A.13) 



-2 sin 9(d - 2 ) exp( — r-^z)da = I m+IV 

a J 0c sm(9) 

The remaining integrals are of the form 
J exp(- Vl - a 2 )f(a)da, which can only be solved in 
very special cases, such as f(a) = a. However, these integrals 
can be evaluated numerically in a fast and reliable way with 
the use of simple q uadrature formulae of low order (see, e.g. 
IPress et alJlll986-200ll) \ 

Finally, for 9 > 9 C the entire conttibution I\ vanishes, since 
in this case all slices are of the type '2', where all photons emit- 
ted from the top surface will leave the system through the side. 
So the total photon intensity emitted by a plasma in the form of 
a thin disk is 



1(9 <9 C ) = h(a c (6j) + h(a c (9)) 
1(9 > B c ) = h(a c = 0). 



(A. 14) 
(A. 15) 



involved here is very similar to the S i contributions, 

p^cosfl x 

S2 — 2 I (l-exp(-a ))dx 

Jo ' sin 9 w(a) cos 9 

+(d - h(a)) sin 9(1 - exp(-T vv(a) / sin#)) (A.ll) 
= 2w(a)cos6' 

+(1 -exp(-T w(fl) /sin6i))- 



Appendix B: Solving the radiation transport 
equation with the exact square box geometry 

The exact solution of the radiation transport equation takes a 
much simpler expression if one uses a box geometry. In fact, 
for this geometry all integrals are of the same type as the region 
called ' 1 ' in the last appendix, which does not involve any kind 
of numerical integrations. To keep the volume of the system 
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constant, we use a = yjnr for the quadratic top and bottom 
surfaces and get 

hoxiO <9 C ) = a(l - exp(T^/ cos 0)) ■ 

, , . sin cos 
(a sin - a cos - 2 ) 

a 

+2adsin0 (B.l) 

hox(0 > 6c) = a{\- exp(r fl / sin0)) • 

. sin0cos0 
(a cos - d sin - 2 ) 

a 

+2a 2 cos0. (B.2) 
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